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1.  Introduction 

The  purpose  of  this  paper  is  to  put  together  in  a  concrete 
form  the  mathematical  theory  of  an  inverse  scattering  problem  which 
was  originally  formulated  in  quantum  mechanical  terms  but  which  is 
related  to  the  problem  of  determining  the  variation  of  electron 
density  in  an  ionized  gas  such  as  the  earth's  ionosphere  by  means  of 
radio  sounding  experiments.   The  theory  described  here  was  developed 

by  H.  E,  Moses  and  the  present  writer  and  was  originally  presented 

[1  2l 
rather  generally  and  abstractly.   '    In  the  present  article  the 

definitions  and  proofs  of  theorems  have  been  revised  in  an  attempt  to 

make  the  analysis  understandable  to  a  reader  who  may  not  be  familiar 

with  the  terminology  of  quantum  mechanics  or  of  the  associated  Hilbert 

space  theory.   While  the  usual  terminology  employed  by  physicists  has 

not  been  discarded  completely  here,  definitions  and  relations  are  given 

in  concrete  terms  whenever  possible,  and  abstract  relations  are  used 

for  convenience  of  notation  rather  than  for  the  purpose  of  obtaining 

great  generality. 

The  physical  problem  at  which  the  present  work  is  aimed  is  to 
find  the  variation  of  electron  density  in  a  weakly  ionized  gas  from 
a  knowledge  of  the  scattering  amplitude  resulting  from  the  incidence  of 
a  plane  electromagnetic  wave.   It  is  assumed  that  the  relevant  properties 
of  the  medi\m  vary  slowly  enough  so  that  a  scalar  wave  describes  the 
scattering  phenomena  with  sufficient  accuracy.   It  is  also  asstimed  that 
losses,  ordinarily  attributed  to  the  collisions  of  the  electrons  with 


heaAAy  particles  may  be  neglected.   The  effect  of  an  external  magnetic 
field  which  would  cause  anisotropy  is  neglected  as  well,  so  that  the 
scattering  is  asstimed  to  be  isotropic  as  well  as  lossless. 

A  steady-state  scalar  field  component  u(x,y,z)  under  these 
conditions  will  satisfy  a  differential  equation 

(A)  Au  +  Jk^  -  V(x,y,z)L  u  =  0, 


having  the  same  form  as  the  time-independent  Schrodinger  equation  of 
quantum  mechanics.  "^   The  quantity  k  is  the  wave  number  and  the  function 
V(x,y,z)  depends  upon  the  variation  of  electron  density  in  the  medi-um; 
in  fact,  V(x,y,z)  is  the  plasma  frequency  divided  by  the  square  of  the 
velocity  of  light  in  vacuum. 

The  scattering  problem  associated  with  such  a  differential 
equation  would  be  to  calculate  the  scattering  amplitude  T(x,y,z,k) 
defined  by  the  asymptotic  form  of  a  solution  u  for  which  the  incident 
wave  is  plane  and  uniform  and  propagates  in  the  direction  of  the  vector 
k.   That  is,  if 

ik'X 

u.     =  e , 

inc. 

where  x  is  the  radius  vector  whose  components  are  x,y,z,  and  at  large 
distances 

ikx 
(B)  u  ~  e^^"^  +  T(x,k)  ^  , 
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where  k  and  x  are  the  magnitudes  of  the  vectors  k  and  x,   the  inverse 
scattering  problem  woiild  he  to  determine  the  function  V(x,y,z)  in  (A), 
given  the  scattering  amplitude  T(x,k),  or  some  part  of  it. 

At  the  present  time  the  inverse  scattering  problem  for  the 
system  described  precisely  by  equation  (A),  in  three  dimensions ,  is 
still  open.   In  this  paper  we  shall  consider,  instead  of  (A),  the 
more  general  integro-differential  equation 

(C)  Z^su  +  k  u  -  /  v(x,x' )u(x' ,k)dx'  =  0, 

where  the  integration  as  indicated  by  the  differential  dx'  is  over 
all  three-dimensional  x'  space.   We  assume  that  the  integral  operator 
whose  kernel  is  v(x,x'  )  is  such  that  the  asymptotic  relation  (B) 
still  holds  for  some  solution  of  (C).   Physically,  the  differential 
equation  (C)  would  govern  electromagnetic  wave  propagation  in  a 
soTiewhat  artificial  plasma  of  a  more  general  type  than  that 
characterized  by  (A).   However,  we  allow  the  kernel  v(x,x'  )  to  be 
a  distribution,  and  in  particular  it  could  have  the  form 

(D)  v(x,x'  )   =  V(x)  5(x-x'  ), 

whereupon  the  integro-differential  equation  (C)  would  be  indentical 
with  the  differential  equation  (A). 

From  the  point  of  view  of  the  inverse  scattering  problem  we 
woiild  have  to  discover  a  method  of  solution  which  would  automatically 
impose  the  condition  (D)  if  we  wished  to  solve  the  original  problem 


associated  with  (A).   Instead,  what  we  actually  do  here  is  to  impose 
conditions  that  guarantee  a  unique  solution  of  the  inverse  problem 
associated  with  (C)  but  which  do  not  guarantee  the  form  (D)  in  general, 
although  the  possibility  is  not  excluded  that  for  some  given  scattering 
operators  the  solution  v(x,x' )  will  have  the  form  (D)  and  hence  will 
apply  to  the  problem  associated  with  (A). 

In  section  2  of  this  article  various  conventions  which  are  used 
throughout  are  listed  for  convenience  and  are  referred  to  specifically 
later  on  whenever  it  seems  appropriate  to  do  so. 

In  section  3  an  integral  transform  theorem  analogous  to  the 
Fourier  integral  theorem  and  involving  arbitrary  solutions  u(x,k)  of 
(C)  is  presented.   This  theorem  requires  the  use  of  a  certain  weight 
fionction  w(k,k' )  needed  to  normalize  the  particular  function  u(x,k) 
which  appears  in  each  instance  of  the  theorem.   We  calculate  w(k,k' ) 
for  the  case  in  which  u(x,k)  is  the  function  satisfying  the  asymptotic 
condition  (B).   The  discussion  should  be  of  some  interest  by  itself 
for  other  more  usual  applications  such  as  the  formation  of  wave  packets 
in  the  medium  determined  by  (C)  (or,  in  particular,  by  (A))  associated 
with  normalized  scattered  plane  waves. 

Section  (h)   deals  with  an  inverse  scattering  problem  for  (C) 
wherein  conditions  sufficient  to  guarantee  uniqueness  of  the  solution 
are  imposed.   It  is  hoped  that  the  same  or  a  similar  formal  setup  can 
be  used  in  the  future  to  attack  the  inverse  scattering  problem  asso- 
ciated directly  with  (A). 
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2.  Notation 

In  this  section  we  list  a  set  of  conventions  for  easy  reference 
later  on.   The  notation  described  here  will  be  used  consistently 
throughout  the  rest  of  this  article. 

a)  An  underlined  letter  such  as  x  will  represent  a  three-dimensional 
vector  (x, ,x^,x_).   The  same  letter  x,  not  underlined,  will  represent 


/  2   2   2 
the  magnitude  V  x,  +Xp+x   of  that  vector. 


2  "3 

b)  Unless  there  are  explicit  limits,  an  integral  will  be  over  the 
entire  domain  of  the  variable  occurring  in  the  differential.   Thus 

OO  00    00    00 

/  f(x)dx  means   /  f(x)dx,  while   /  f(x)djc  means   /   /   /  f  ( x^ ,  x^ ,  x^  )dx^' 

O  -00   -OO   -00 

Miatiple  integrals  will  be  designated  by  the  nature  of  the  differentials 
which  occur  in  them,  and  a  single  integral  sign  will  always  be  used. 
For  example,  /  f(xjZ)<3^^^  is  ^  six-dimensional  integral  over  the  entire 
domain  of  the  vectors  x  and  ^. 

c )  We  shall  be  concerned  with  linear  operators  in  the  form  of 
integral  operators  over  a  suitable  space  of  functions.   The  kernel  of 
such  an  operator  will  be  denoted  by  a  lower  case  letter.   When  we 
wish  to  speak  of  the  operator  in  the  abstract,  the  corresponding 
capital  letter  will  be  used.   Thus,  Uf  means  the  abstract  operator  U 
acting  on  the  abstract  linear  element  f ,  and  the  concrete  representer 
of  U  will  be  /  u(x,x'  )f(x'  )dx' . 


dXgdx  . 


d)  The  kernels  allowed  for  integral  operators  can  be  distri- 
butions as  well  as  smooth  fionctions.   Thus,  in  particular,  the 
identity  operator  I  is  represented  by  a  kernel  which  is  a  delta  function; 


'I 


If  =  /  S(x  -x'  )f(x'  )djc'  =  f(x)  =   f. 

The  Laplacian  operator  A  can  be  represented  by  an  integral  operator 
whose  kernel  is  given  by  A&(x-x'  )  =  A'6(x-  x'  ),   wherein  the  primed 
differential  operator  is  taken  with  respect  to  the  primed  variables, 
and  the  unprimed  operator  with  respect  to  the  unprimed  variables.   For 
example,  we  have 

Af  =  Af(x)  =   /  Ja6(x -x'  )l  f(x'  )djc'. 


e)  ¥e  shall  use  the  notation  U*  for  the  Hermitian  adjoint  of  U, 
but  u*(x,x'  )  will  mean  the  complex  conjugate  of  the  function  u(x,x'  ). 
Hence  the  following  are  equivalent: 

U*f  =   /  u*(x',x)f(x' )dx'  . 


=   /  u*(x' ,x 


f )  Instead  of  functions  of  x  the  Fourier  transforms  of  the 
functions  will  sometimes  be  used  to  represent  the  same  linear  space. 
The  new  functions,  the  transforms,  will  depend  on  vectors  k  instead  of 
X  and  will  be  indicated  by  a  circumflex.   Thus 


f=f(x),   f  =  f(k)  =— ^-^  y^e^^'^  f(x)d^. 


(2«)- 


In  the  case  of  an  integral  operator  kernel  we  shall  also  write 


Because  of  the  identity 


(2rt) 
we  have 


_J^      f    3iiS-(x-x')  ^  =  S(x-x') 
(2rt)3  J 

ru(x,x' )f(x' )di'  =  /  u(x,k)  f(k)dk. 


g)  The  Laplaclan  operator  A  which  acts  on  the  functions  of  f (x) 
is  represented  in  k-space,  the  space  of  the  Fourier  transforms  f(k), 

by  the  multiplicative  operator  -k  (cf.  2. a),  that  is,  by  the  integral 

2 
operator  whose  kernel  is  -k  S(k-  k' ).   Thus, 


Af(x)  =  -  ^  Tk^e'^^'^  f(k)dS  =  -   ^  .<     fy^^Hh-  k'  )e"^-^  f(k'  )dk'dJs. 

(2«)3/2j  (2::)3/2J 


Hence  we  can  write 


A  =  -k^  =  -k^S(k- k'  ) 


k-space.   Because  of  the  relation 
f(x)&(x)  =  f(0)6(x) 


any  operator  whose  kernel  has  the  form 


r  =   r(k  -k'  )  &(k^-  k'^) 


commutes  with  the  Laplacian  A.   That  iSj 


AT  =  TA, 


for  we  have 


-k^r(k-is')&(k^-  k'^)  =  -k'^r(k- k')&(k^-  k'^). 


3.  Eigenfimctions 

Solutions  u  (x^k)  °f  'tke  differential  equation 

(1)       -A  u  (x,k)  =  k^u^(x,k) 

will  be  called  eigenfunctions  of  the  operator  -A,   Each  solution 
u  (x,k)  will  ^e  regarded  as  the  kernel  of  an  integral  operator  which 
is  a  representer  of  some  abstract  operator  U  .   This  integral  operator 
acts  on  the  space  of  functions  f(k)  which  are  transformed  by  it  into 
functions  f(x)  belonging  to  the  x-representation.   Thus,  we  shall  have 
relations  of  the  type 

cp(x)  =   /  u^(x,k)  f(k)dk. 

A  special  case  of  this  is  the  Fourier  transform  for  which 
/   1  N      1      -ik-x 

(cf.  2.f). 
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Similarly,  we  shall  be  concerned  with  operators  U  and 
corresponding  eigenfionctions  u(x,k)  which  are  solutions  of  the 
integro-differential  equation 

(2)       -Au(x,k)  +   /  v(x,x' )u(x' ,JS)^'  =  ^   u(x,k). 

In  operator  form  (2)  is 

-AU  +  VU  =  k^U. 

Here  again  u(x,k)  is  regarded  as  the  kernel  of  an  integral  operator 
which  transforms  a  function  of  k  into  a  function  of  x-   In  (2)  the 
operator  V  whose  kernel  is  v(xjX' )  "^^y  in  special  cases  be  a  multi- 
plicative operator;  i.e.,  v(x,x' )  may  have  the  form  cp(x)&(x-x') 
(cf.  2.d).   Then  (2)  will  have  the  usual  form  of  the  Schrodinger  equation. 

It  is  well  known  that  the  operator  -A  is  self -ad joint  and  has 
a  continuous  spectrum  consisting  of  zero  and  all  positive  real  numbers 
'when  acting  on  a  space,  e.g.,  of  square  integrable  functions  f(x).  We 
shall  assume  that  this  is  true  of  -A  and  in  addition  of  -A  +  V  without 
worrying  about  particular  conditions  on  the  "perturbation"  V  necessary 
to  guarantee  this  property  or,  for  that  matter,  about  the  precise 
nature  of  the  space  of  functions  f('x). 

We  have  given  an  operator  interpretation  of  the  eigenfunctions 
of  (l)  and  (2)  as  transformations  from  k  space  to  x  space.   The  more 
concrete  interpretation  is  also  useful,  especially  in  applications,  in 
which  we  regard  them  as  providing  integral  transform  theorems  analogous 


to  the  Fourier  integral  theorem.   We  can  then  use  them,  for  example, 
to  describe  the  transient  response  to  the  system  described  by  (l)  or 
(2).   Then,  the  Fourier  integral  theorem  provides 


(3)  .         f(x)  -  — ^  /  e"^-"-  f(k)dJv, 


Ik'X  C;, 

77^   /   e 
(2rt)- 


(M  f(k)  =  ^   /   e^^-^  f(x)d^. 


-ik-x 


I^^ 


then  (3)  and  [h)   imply  the  operator  relation 

(5)  U  U*  =  U*U  =  I, 
^    '  o  o    o  o 

where  I  is  the  identity.   The  relations  (5)  mean  that  U  is  unitary. 

Similarly  if  we  define  the  particular  eigenfunction  u  (x,k) 
of  (2)  as  the  unique  solution  of  the  integral  equation 

(6)  u^(x,k)  =  ^-i^  e^^-^  -  ^  /^^^^f^  v(x',x")u^(xMs)d.'d.'V 
then  we  can  prove  an  analogous  transform  theorem  in  terms  of  u  (x,k) 


The  transient  response  would  be  given  by  multiplying  a  function  f(k; 
by  e^^u(x,k),  where  w  =ck,  and  integrating  over  w. 

Tlie  solution  of  the  integral  equation  (6)  is  unique  because  of  our 
assumptions  about  the  operator  -A  +  V. 
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and  u^(x,k).   We  can  prove: 

if 

(7)  f(x)  =  j^u^(x,k)9(k)dk, 

then 

(8)  cp(k)  =  r  u*(x,k)f(x)(bc. 

Relations  (7)  and  (8)  imply  the  relation 


(9) 

for  the  operator  U  corresponding  to  u_^(x,k).   Thus,  U_^  is  imitary. 

Similar  resiilts  hold  for  the  eigenfimction  u_(x,k)  and  the 
corresponding  operator  U  ,  where  u  (xjk)  is  defined  as  the  solution  of 
the  integral  equation 

(10)      u  (x,k)  =   ^  ^/p  e^-"-  -  -^  f =^-=^-  v(x',x")u  (x",k)d^'d^" 

(2rt)-^/  ^    |x-  x'  1 


We  shall  prove  the  assertion  we  have  made  concerning  the  eigen- 
functions  u^(x,k)  by  showing  that,  for  example, 


(11)      N*(x,k)u_j_(x,ls'  )djc  =  6(k  -k'  ). 


Since  we  have  assumed  that  the  spectra  of  -A  and  -A  +  V  are  the  same 
the  self-adjointness  of  these  operators  will  imply  that  U  has  a  right 
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inverse  as  well  as  a  left  one.   Thus,  "because  of  (11)  we  shall  have 

uX  =  ^<  =  I. 

The  arg-uments  required  to  prove  a  similar  theorem  in  terms  of  U  are 
completely  analogous,  and  we  shall  have  finally 


(12)  U,U,  =  U.U,  =  I, 


(13) 


/  u+(2S.k)u+(x,k'  )dx  =  &(k  -k'  )> 
/  u+(x>k)u+(x'  ,k)dk  =  &(x-  x'  ). 


Proof  of  (11):   Before  beginning  the  argument  it  will  be  usef\il 
to  define 

4f(xjk)  -      /  v(x,x' )u  (x' ,k)djc' . 
We  can  now  write  in  place  of  (6) 

We  shall  also  need  the  well-known  identity 

±ik|2i-x'  I       ,.      ^     p   i£-(x-x') 
/•nc\        1  e    1-^  -  I       iim     1    /  e  -^  ' , 


-  13 


which  is  proved  in  appendix  I.   The  relation  (15)  has  the  meaning  of 
a  distribution;  i.e.,  the  limit  is  to  be  taken  after  an  integration 
over  the  variable  x ' : 


We  shall  also  require  the  identity 


^-■-^^    (a±ie)(b±ie)   (a+b)±2ie  ]  a  ±ie  "^  b  ±ie 


which  follows  from  a  simple  straightforward  algebraic  manipulation. 
In  addition  we  shall  need 

r    -ik-x +  ik'  Ix-x'  I  -.  .  r    ^~^^'^' 


and 


(18) 


P  ik' -2;  -  ikl^i-x' I  -,  .    n      ik'-x' 

t      [^ —=-        t"(xSk)d^'dx=  ^^"   r  g     .     A2cMs)d^- 


Relations  (17)  and  (l8)  follow  easily  from  (15).   Thus,  for  example 

-h    /^"^-■-^'''''-"-'^(2c',k')dx'dx 

lim  1  r    iWx-x')-  ik-x 

-     e^o       77~3     /  ; 2 t(x',k' Mx'dx  dp 

[2n)      d        k'^-  k    +  it 

^    k'     -  p    +  ie 

-,  .  p        -ik'x' 

"e^O       /  ~;2 2 t(x',k')d^'. 

^      J    k      -k   +  ie 
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Finally,  we  shall  require  the  identity 


(19)      ru*(i,k)t^(ir,k'  )di  =  J u{x,K'  Hl(x,h 


)di[. 


This  relation  follows  from  oijr  assumption  that  the  operator  V  is 
Hermitian,  i.e., 

v(x,i)  =  V  (i[,x). 
We  have,  in  fact, 

/  u*(:i:,k)t(z.k'  )<iii  =  /  u*(y:,k)v(z,2i)u(i,k'  )dxd^ 

=  /  u(x,k'  )v*(x,ii)u*(2,k)djcdx  =  /  u(x,k'  H*(x,k)dx. 

To  complete  the  proof  of  (11 )  we  apply  the  integral  equation  (ih), 
using  the  right  side  in  place  of  u  (2c,k)  and  its  complex  conjugate  in 
place  of  u*(x,k),  to  foi*m 


/  u*(x,k)u^(x,k'  )djc  =  &(k  -k'  ) 


(20)  p  -ik-x  +  ik'  |x-x'  [ 


375  J — ;~7' — v<2''i''^'^ 


(l|jr)(2rt)-^'    ^       |x  -X 


£i.    I 


P  ik' -x  -  ik|x-x'  I 

^^—j-      /^ t,(x',k)'ix'd2c 

(l^rt)(2rt)^/"^  ^      |x  -X'  I 


-,  p  ik'  |x-z'  I  -ik|x-x'  I 

+  — i-  /^ <(x',k)t,(:z',k')djc'di:'dx. 

(hnfJ  |x-:^'|Ix-x'| 
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We  now  consider  the  last  term  on  the  right  of  (20)  and  apply 
(15)  to  the  functions  which  are  factors  of  \(f  and  \|/  in  its  integrand. 


1    Pg-ik|x-x'|  +ik'|x-i[' I   ^ 
{k^f    J  |x-x'||x-i'| 

Tim    1  r     pi£"(2c-x')+i^-(x-i[' )   ^ 

=  e  ^0  ~^     /  ~T"^ rrri ? <(^'  'is)n(:i'  .k'  )d2c'di'dxd£d£ 

^  -"^  (2rt)^  ^  (k^-p^-  ie)(k  2_q2^.^^   +       + 


Tim     1      r   p-^-2^  "  ^^'^ 
-  -07773  J  J  2  .  ,,„2  2  .  ,  <(2Sk)t,(x',]i')5(E+a)d£d^d.'dx' 
(2rt)-^  ^  (k  -V   -iejfk'  -q  +ie) 


lim    1     /•    ^^-i^'-^'^ 


n        P       „!£•  {X     -X-  J 

n  -^  /  —^ ^-^ K^^'  ^i^H.  (x'  ,k'  )d£d^'dY' 

°  (2ny    J    (k%2-ie)(k'2-p2+i€)  ^       "^ 


^0-^  T-tV -T-i <(x'.is)n(zMs')d£d^'d:^'. 

-*0  rp^^3  J  ,  2  2    w  ,2  2    .   +       + 


^  """^  (2jr)-^  ^  (p"-k"+ie)(k'"-p"+ie) 


If  we  apply  the  identity  (I6)  this  expression  hecones 


lim 
e  ^  0 


-lim 


my    J    k^-k^+2i£    I-'      k^-p^-ie  ^  J      "" 

^     f—^h 1    r^!?     2     ~       t,(x',k')<iy'd£U:!(x',k)djc' 

!rt)3    -'    k'2-k  +2ie    y    k'^-p  +ie  "^  J      "^ 


(2rt)-^    ^    k     -k  +2ie    L^    k     -p  +iG 
If  we  now  \ise   (17)  and   (I8)  it  becomes 
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-lira 
e  -^  0 


We  now  substitute  from  (3U)  for  the  bracketed  expressions  and 
obtain 


-lira  n  '        --^^-^ 


^  -0     j,.2_j^.^2ie    ^       (2« 


+li 
e 


1  r      e'  -^ 

k2+2ie    ^       (2rt)3/2     + 

^To-rh ru:(z',k)V^',k')dz' 

"^k*    -k  +2ie         ^ 

-,  p        ik'  -x' 

^°      k'2-k2+2ie     ^        (2:t)3/2         + 

^  k'^'-k  +2ie   ^ 


-lim 


Wow  by  (19)  the  second  and  fourth  terms  in  this  expression  cancel 
each  other.   Hence  we  have  finally  for  the  last  term  on  the  right  side  of 
(20)  the  expression 


-lim     ^        r.-ik-^' 


k'2-k^+2ie  ^  (2Tr)3/2  + 


^-°— —  -     (2.)- 


lt  •       T       P  ik'  -x' 
^  ^"  k'"^-k  +2i6  ^  (2rt)^/^ 
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If  we  make  use  of  (1?)  and  (l8)  we  see  that  the  two  terms  in  this 
expression  cancel  the  second  and  third  terms  on  the  right  of  (20), 
and  the  relation  (11 )  which  we  set  out  to  prove  now  follows. 

We  can  define  more  general  eigenf unctions  of  -A  and,  corres- 
pondingly, more  general  eigenfunctions  of  -A  +V  can  be  defined  by- 
means  of  integral  equations  similar  to  (10).   For  example,  we  can 
define  the  eigenfunction 

(21)  vi^(x,k)=   ^   r  e^^-^m_^(k',k)dJc' 


of  -A.   The  function  u   (x,k)  will  be  an  eigenfunction  of  -A  if  we 
assume  that  m,(k',k)  has  the  form 


(21a)         m^(k',k)  =  u^(lS',k)  S(k'2-k2), 

according  to  the  remarks  in  (2.g). 

We  can  then  define  a  corresponding  eigenfunction  u(x,k)  of 
-A+V  as  the  solution  of  the  integral  equation 


ik|x-x' I 
(22)     u(x,k)  =  u^_^(x,k)  -  fz      I   "^  i    ^,  i   v(x' ,x"  )u(x",k  )dx'djc 


o+' 


W  J 


hn    J      |x-x'| 


If  we  replace  k  by  k'  in  (10),  multiply  through  by  m  (k,k'  ) 
and  integrate  over  k',  we  see  fron  the  "uniqueness  of  the  solution 
of  (22)  and  of  (10)  that 
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(23)      u(x,k)   =   /  u^(x,k' )m^(k',k)dk'. 
Expressed  in  terms  of  abstract  operators  (23)  is 

(2U)  U  =  U^M_^, 

wherein  U  is  the  operator  corresponding  to  the  kernel  u(x,k)  and  M_^ 
is  the  operator  corresponding  to  the  kernel  m  (kjk'  )• 

We  woTjld  also  like  to  define  an  operator  M_  such  that 

(25)  U  =  U_M_. 

We  can  do  this  by  showing  that  the  eigenfimction  u(x,k)  also  satisfies 
an  integral  equation 

p      -ik|2c.-x'  I 

(26)  u(x,k)  =  u^_ix,k)--^      /  ^  Ix-x'l   ^(j?'^2S")^(x'M0d-x'dx", 

where  u  _(x^k)  is  an  eigenfunction  of  -A  given  by 

(27)  u   (x,k)  =     ^,/p   r  e^^'"^m  (k',k)dk'. 

°"         (2n)3/2  J 

Here,  again,  m_(k',k)  must  have  the  form 

(27a)    m_(k',k)  =  ^_(k' ,k)S(k' 2  -  k^ ). 

We  can  obtain  a  relation  for  m_(k',k)  by  setting  the  expression 
(26)  for  u(x,k)  equal  to  the  expression  (22).   We  have  then 
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^    p    ik|x-x' I    -ik|x-x' | 
(28)    u„.(x.l< )  .  u„^(x,!s)  -  i^  J     p^^   -  °U-;.j 

X  v(x' ,x" )u(x",k)dx'dx". 
-ik' -x 

We  now  multiply  (28)  by  ryr  and  integrate  with  respect  to  x. 

(2n)3/2 


According  to   (27)  and   (23)  the  resist  is 


(29) 


,         p      -ik'-x      I    ik|x-x' I  -ik|x-x'|] 

m.(k',t)  .  ^,{K',K)--,   J     ^^     |ni^  -   '    U-Vl] 


^  v(x' ,x"  )"(x'Ms)d^'^i"<^- 


It  is  clear  from  (28)  and  the  properties  of  the  Green's  functions 

±ik(x-x'  ) 
— ] n —  that  V.     (x,k)  is  an  eigenfunction  of  -A.   Hence,  the 

form  (27a)  of  m  (k',k)  we  have  assumed  must  be  correct. 

We  can  prove  this  more  directly  from  (29)  by  making  use  of  (17) 
and  (l8).   Thus,  in  place  of  (29)  we  have 

(30)  m  (k',js)  =m  (k',k)+^^"^^   ^^/g   /  {% ^ --^ W(x',k)cbc' 

""        ^^°(2jr)3/^^   [k^-k'+ie   k^-k'^-iej 

■ik'  -x' 


m.(k',k)— ^  r:^o^''f  X'2^2,.2-^^^'-^^^' 


1      lim 
(2^     ^  -"^  d       (k"-k'")%e' 


/-ik' • x' 
^     -     -     t(x',k)d2c'6(k'2-k^), 


{2«)3/^ 


where  we  have  used  the  Identity 


X    +  e 


The  form  (27a)  for  m_(k',k)  follows  from  (30 )  and  the  corresponding 
form  (21a)  for  m_^(k',k). 

If  either  one  of  the  operators  U  or  M  is  known  to  have  an  inverse 
then  (2U) 'implies  the  existence  of  the  inverse  of  the  other.   Then 
from  (25)  the  existence  of  the  inverse  of  M  follows.   If  in  (2k)   U 
is  set  equal  to  U  then 


In  this  case  (30)  becomes 

fe-^^'-^'  .2   2 

(32)  m  (k',k)  =  &(k'-  k)- 27Ci  /  ^ TTT"  ^^(xMs)d^'S(k  -  k  )  =  S, 

^  (2n)3/2   + 

The  operator  S  defined,  by  (32)  is  called  the  scattering  operator  (from 
the  definition  employed  in  quantum  mechanics)^  and  from  what  we  have 
just  observed  it  has  the  form 

(32a)   S  =   s(k',k)  =  a(k',k)5(k'2-k2). 

The  forms  (21a),  (27a)  and  (32a)  imply  in  each  case,  we  recall  (cf.  2.g)j 
that  the  operator  commutes  with  the  Laplacian  A. 
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The  relationship  between  the  kernel  of  S  given  in  (32)  and 
the  asymptotic  value  of  u    (x^ls)^  which  may  "be  interpreted  as  the 
scattering  amplitude  for  an  incident  plane  wave,  can  be  seen  at  once 
from  a  consideration  of  (6)  in  the  limit  of  large  x.   We  observe  that 


2  o    .  ,   r2 


I X  -  X '  I  =  y/x~ -   2x •  X '  +  X '   ~  X 


where  x  is  a  unit  vector  in  the  direction  of  x.   Then  from  (6), 
-o  - 

-,  /^      -ikx  -x'  ., 

ik-x    /^  v3/2  n  -o  -  ikx 

(33)    u^U,K}  ~  ^^  -  ^^ft-    J     7^;;^  v(x.,x")uJ_xM.)d..c^"_. 

ikx 
The  factor  of  .§ in  the  second  term  on  the  right  of  (33)  is  the 

f27r)3/2 
scattering  amplitude  -  ^"-j^ t_^(k',k),  where 


/-ik' -x' 
TT^  v(x',x")u,(x",k)ci^'d2c"  > 
(27t)3/2  + 


and  it  is  assumed  that  k'=k  and  k'  has  the  direction  of  x-   Then  we 
can  write  for  the  kernel  s(k',k)  of  S  the  relation 


(35)   s(k',k)  =  5(k'-  k)  -2jtit^(k',k)S(k'^-  k^). 


We  have  from  (2^),  (32)  and  the  manner  of  deriving  the  relation 
for  S  that 

(36)  U^  =  U  S, 
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and  hence  that 

(37)  S  =  uX. 

We  have  used  the  fact  that  U_  is  imitary  in  deriving  (37),  which  in 
turn  implies  that  S  is  unitary: 

ss*  =  (u*u^)(uX)*  =  U*(U_^U*)U_=  U*U_  =  I. 

We  can  conclude  from  this  that  S  is  unitary  because  we  have  assumed 
that  U  and  U_  have  left  and  right  inverses. 

From  (2l+)  and  (25)  we  have 

U  M  =  U  M_  . 

Thus, 

M_  =  U^^M^  =  SM^  , 

from  which  we  get 


(38) 


We  can  now  state  an  integral  transform  theorem  in  terros  of  an 
arbitrary  eigenfunction  u(x,k)  of  the  operator  -A  +V. 

Given  an  arbitrary  function   f (x)  of  the  type  we  consider  as 
representing  the  abstract  linear  space  (e.g.,  square  integrable  functions) 
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we  can  write 


(39)  f(x)     =    f   u(x,k)cp(k)dk, 

where 

k)     =     /  f(x)  I      /  w(k,k'  )u*(x,k'  )dJc'  >  dx 


(UO)  cpQ 


and  the  ftinction 

/  w(k,k'  )u(x,k'  )dk' 


is  an  eigenf unction  of  -A  +V.   This  can  be  guaranteed  by  reqiiiring 
that  w(k,k' )  have  the  form 


(U)         w(k,k')  =  uj(k,k' )5(k2-  k'^). 


It  will  also  turn  out  that  the  operator  W,  which  we  shall  refer  to  as 
the  spectral  weight  operator,  corresponding  to  w(k,k' ),  is  Hermitian 
and  positive  definite. 

We  can  define  W  abstractly  by 
(1^2)  W  =  m;^^*"^  ; 

then  from  {2k)   and  (25)  we  obtain 
(U3)  W  =  m"^*"^  . 
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The  fact  that  W  as  defined  by  (U2)  or  (^3)  is  Hermitian  and  positive 
definite  is  obvious. 

The  relations  (2l|)  and  (25)  also  imply 

(UU)  WU*U  =  I 

and 
(U5)  UWU*  =  I. 

Relations  (hk)   and  (ij-5)  are  equivalent  to  the  transform  theorem  relations 
(39)  and  (UO). 

From  (k2)   or  ('^3)  we  see  that  the  spectral  weight  operator  W 
needed  for  the  integral  transform  theorem  corresponding  to  a  given 
eigenfunction  u(x,k)  can  be  determined  if  we  know  the  integral  equation 
of  the  form  (22)  which  u(x,k)  satisfies.   We  can  also  obtain  the 
spectral  weight  operator  W  in  terms  of  the  asymptotic  amplitude  of 
u(x,k)  as  X  becomes  large.   The  scattering  operator  S  plays  a  signi- 
ficant role  in  this  relationship,  and  we  can  therefore  obtain  a  direct 
connection  between  S  and  W. 

Consider,  for  example,  an  eigenfunction 
ik6)  u(x,k)  =   ru^(x,k'  )m^(k',k)dk' 

defined  to  have  the  asymptotic  behavior 


,   ,  V    ik-x 
u(x,kj  ~  e  


as  x  becomes  large  in  a  "forward"  direction,  that  is,  for  the  polar  angle 


25 


e  of  X  in  0  <  9  <  |.   We  can  write  instead  of  (h6) ,   vith  m_^   (k',k) 
as  the  kernel  corresponding  to  M  , 


(hi 


)  ru(x,k' )ra;^-^(k',is)dis'  =  u_^(x,k). 


We  now  consider  the  asymptotic  form  of  (U7)  for  large  x,  using  (33) 
and  (3^)  and  (A.l)  of  Appendix  II.   We  have,  assiming  that  m   (k'jk) 


has  the  form 


r^(k',k)  =  n;^(k',k)s(k'^-  k^). 


(i+8)         fvi(x,k<)ml^(k',^)dk'~      re^-'"V-^(k',k)dJs' 


o   o   o   ^ 


2jr     jt 


^  J   dcp'  TdG'  rdk'M^^(k',k) 


^       o      o     o 


ikoc  -ikx 

X  ^-6(k'-k)6(©'-e)6(cp'-cp) +&(k'-k)S(jt-e-0'  )5(9+jt-cp'  )^ 


(2Tt) 


37ii"^^^^o'^)V-"^^^-^o'?^)^''" 


where  k  is  the  vector  -  x.   In  the  calculation  of  (i+8)  we  have  used  the 
identity 

S(k'^-k2)  =  5([k'+k][k'  -k])  =^  6(k'  -k). 
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We  also  have  from  (33) 


ikx  -ikx 

(1^9)  u^(x,k)  ~— ^^^-6(e-e^)6(cp-cp^)^  +  S(n-0-G^)&(cp±«-cp^)  ^-^ 

^ t  (k  ,k)^—  . 


By  combining  (^7),  (^8)  and  {h9)   we  obtain 


^:^(k  ,is)  --^ 5(9-9  )S(cp-cp^)-2jtit(iS_,k), 

(50)    ^   °      k^sin  9       °       °       ^ 

1^+  (-k^.k)  =  -:^ 6(rt-  9-0^)6(cp±jr-cp^). 

k  sin  9 


Finally,  from  (35),  (50)  and  the  identity 


5(1,,  _  1,)  ^  5(k'-k)5(e'-9)5(cp'-9; 
k  sin  0 


we  have 

"i+^(ko'is)  =  6(k^-  k)-2nit^(]i^,is)6(k^-k^)  -  s(k^,k),   0  <  9^  <  | 
(51) 

m^"^(-kQ,k)  =  S(kQ+k)  or  ^^ [^^,\)   =  HK^-   k),  |  <  9  <  jt. 


(52)      m;l(k^,k)  =  s(k^,k),  (I  -9J   .  6(k^-  k),  (e^-l) 
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Then  from  (^1^2)  we  obtain  the  kernel  corresponding  to  W  in  the  form 
w(kQ,k)  =   rm;j;-^(kQ,k' )ra*"^(k,k' )cik' 

=  /s(ko'k'  )s*(k,k'  )dJs'Ti(|  -eo)Ti|-  9)+  sas^,k)Ti(|-  e2)Ti(0-g 

(53) 

-s^(is,kjn(f-9).(e^-f).5(k^-k),(;e^-|) 

=  5(iSo-is)-^s(is^,kh(|-eJ,(9-|).s*(k,is,h(|-e),(9^-|), 

wherein  we  have  used  the  unitary  property  of  the  scattering  operator  S. 

k.    The  inverse  problem 

In  the  inverse  problem  our   goal  is  to  obtain  the  operator  V 
of  (2),  that  is,  the  kernel  v(x,x' ),  from  a  knowledge  of  the  scattering 
operator  S,  or  actually  from  a  knowledge  of  some  part  of  S.   It  will 
be  convenient  to  rephrase  the  problem  in  a  different  form  in  terms  of 
the  spectral  weight  operator  W.   We  can  do  this  by  making  use  of  (53)j 
and  we  observe  that  w(]i  ,li)  corresponding  to  W  requires  only  part  of  S 
for  its  complete  specification. 

In  the  new  form  the  inverse  problem  is  to  find  V  from  a  knowledge 
of  W.   For  this  purpose  it  will  be  convenient  to  rewrite  (kk)   and  (U5) 
in  the  form 

(5i|)  WU*  =  U 

(55)  U^  =  U"^  . 
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If  a  solution  U  of  (5U)  and  (55)  can  be  foimd  then  V  will  be 
given  by 


'.%) 


u(x,x' )  =   /  u(x,k)  ~ 773  dJc  =   U 

J  (2rt)3/2 


where  here  u(xjx' )  is  the  kernel  of  the  integral  operator  representation 
of  U  on  the  space  of  functions  of  x.   By  (55)  the  result  will  be 

-  /  u(x,x'  )  A'  — — —rr   djc'  =  -k   /  u(x,x'  )  ~ — TT3   djc'  =  -k  u(x,k; 
J  (2:t)3/2  J  (2rt)3/2 

That  is,  we  shall  have  the  expression  equivalent  to  (56): 

-Au  (xjk)+  /  v(x,x'  )u(x',k)(3^'  =  k  u(x,k)  , 

where  v(x,x'  )  is  the  kernel  corresponding  to  the  operator  V  given  by  (56). 

Prom  (^h)   and  (56)  it  follows,  since  A  is  Hermitian,  that  if  V 
is  Hermitian,  i.e., 

(57)  V  -  V*  , 

then 

U  A  WU*  =  UW  A  U*  , 
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where  we  have  used  the  fact  that  W  is  Hermitian.   Since,  according 
"to  (55),  U  has  an  inverse  it  follows  that 

(58)  AW  =  WA  , 

which  means  that  the  kernel  corresponding  to  W  in  the  k  representation 
has  the  form 


(59)  w(k,k')6(k2-  k'^) 


as  we  have  ass\jmed.   Relations  (58)  and  (59)  are,  in  fact,  necessary 
and  sufficient  for  (57)  to  hold. 

Obviously  the  solution  of  (5^),  (55)  for  U  is  not  unique  for 
a  given  W.   As  we  have  just  observed  when  W  has  the  form  (59),  which 
we  always  assume,  then  any  V  corresponding  to  a  solution  U  will  auto- 
matically be  Hermitian.   In  order  to  guarantee  that  V  be  truly 
sejLf-adjoint  and  that  the  correct  scattering  operator  be  reproduced 
we  must  introduce  other  conditions.   The  general  question  of  the 
uniqueness  and  existence  of  V  under  any  general  condition  which  leads 
to  the  correct  scattering  operator  is  still  open,  and  we  shall  not 
attempt  to  settle  it  here.   Instead  we  shall  consider  a  peirticiilar 
condition  which  leads  to  a  unique  V  such  that  the  part  of  the  scattering 
operator  needed  to  define  W,  as  in  (53),  will  be  reproduced. 

The  condition  we  shall  prescribe  is  a  restriction  on  the  form  of 
the  kernel  u(x,x'  )•   We  write,  first, 

(60)  U  =  I  +  K  , 
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and  then  assume  that  K  is  represented  by  a  kernel  of  the  form 

(61)  K  =  q(x,x'  )ti(x'  -  x)  . 

Then 

u(xjX'  )   =  5(x  -  x'  )  +  q.(xjx'  )n(x'  -  x) 

and 

ik-x      p     rs     p  ik-x' 

(62)  u(x,k)  =  -^ — 77^+   /   /   /  q(x,x'  )  -^ — ry-  x' sin  G'de'drp'dx' 

(23t)3/^    J     J     J  (27t)3/2 

^         '  X   o   o  ^ 


From  ( 62 )  we  have  formally 


ik-x 


(63)  u(x,k)  ~     ^/g 


for  large  x  as  required. 

Now  from  (6l)  we  see  that 

{Gh)  K*  =  q(x,x'  )  =  q.*(x' :.x)ti(x  -x'  ). 

We  shall  prove  that  the  assumption  involved  in  (60),  (6l)  and 
{6k)   leads  to  the  following  theorem:* 

For  any  U,  U  of  the  form  given  in  (6o),  (6l)  and  satisfying 
the  equation  (5U), 


The  proof  of  this  theorem  was  furnished  by  Professors  J.  B.  Keller 
and  B.  Friedman. 
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we  have  the  relation 

(65)  ^o  =  ^  ' 

i.e.,  if  U  has  an  inverse  the  equation  (55) 


U  =  U"^ 
o 


will  be  satisfied  automatically. 


PROOF:  We  have 


(66)  UU  =  UWU*  =  U  U^ 

o         o 


since  by  (U2)  w  is  Hermitian.   On  setting 

(67)  U  =  I  +  K  , 

o       o 

where 


(68)  K^  =  q.^{x,x'  hix^  -   x) 


we  have  from  (66) 

(69)  K  +  K^  +  KK^  =  K*  +  K^  +  K^K^ 


32 


In  terms  of  q  and  q  this  is 


q(x,x'  )il(x'-x)  +q  (x,x'  )ti(x'-x)  +  /  q(x,x"  )Ti(x"-x)q^(x",x'  )t](x'-x"  )dx" 

=  ^*(x'  ,x)ti(x-x'  )+  q.*{x',x)T\(x-x'  )+  /  q^(x'  ,x"  )ti(x"-x'  )q*(x",x  )t](x-x"  )dx' 


Each  term  on  the  left  vanishes  for  x'  <  x,  and  each  term  on  the  right 
vanishes  for  x  >  x'.  It  follows  that  both  sides  are  zero  identically. 
Thus  J  we  have  proved  that 


as  desired. 

We  now  observe  that  if  we  set 

w  =  n  +  I, 

then  equation  (5^)  becomes 


m*  +  n  +  K*  =  K 
o 


(70)  Kfi  +  n  +  K  =  K* 


In  terms  of  the  integral  operator  kernels  we  obtain  from  (70),  if  we 

use  the  fact  that  K  and  K  have  the  form  given  by  (6l), 
o 
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(71)    /    /   /  q(x,x")'j(x",x' )x"  sin  e"dx"de"dcp"  + w(x,x'  )+ q(x,x'  )=  0,x'>  X. 

O    O    X 

For  each  fixed  x  equation  (71)  is  an  integral  equation  for  q(x,x' ) 
as  a  function  of  x'  in  the  range  x'  >  x.   If  equation  (71)  has  a  well- 
behaved  solution  for  q(x,x' )  in  the  range  x'  >  x  we  can  form  the  kernel 

u(x^x'  )  =  S(x  -  x'  )  +  q(xjx'  )^(x'  -  x). 

We  can  prove  that  (71)  has  a  unique  solution  by  making  use  of 
the  positive  definite  character  of  the  spectral  weight  operator  W.   It 
follows  then  from  (63)  that  the  solution  we  have  given  reproduces  the 
part  of  the  scattering  operator  required  to  form  the  operator  W. 

We  can  show  that  the  solution  (assumed  to  be  bounded)  q(x,x' ) 
of  (71)  is  unique  by  a  standard  argument  which  goes  as  follows. 
Suppose  q-,  (x^x'  )  is  a  second  solution  of  (71)  and  that  q  and  q^  are 
bounded.   Then  the  function 

r(x,x' )  =  q(x,x' )  -  q-L^^'-'  ^ 


satisfies 

2«   rt 


/   /   /  r(x,x")'^(x"jx'  )x"  sin  9"dx"de"dcp" +r(x,x'  )  =  0,x'>x. 


o   o   x 


3h 


Whence 

/  T)(x"-  x)r(x5x"  )^'^(x"^x' )  +S(2S"  -  x' )}d.x"  =  0. 

This  has  the  operator  form 

rw  =  0. 

We  mtiltiply  (operator  multiplication)  on  the  right  by  P  and  apply 
to  an  arbitrary  vector  cp: 

T\F*(p   =  0. 
The  inner  product  of  cp  with  this  is 

iq>,rw*cp)  =  (r*cp,vr*cp)  =  o. 

Since  W  is  positive  definite  the  vanishing  of  the  quadratic  form  implies 

r*cp  =  0. 
Since  cp  is  arbitrary  then 

r*  =  0, 

and  hence 

q.(x>x' )  =  q-]_(x^x' )  > 

as  required. 

Finally,  we  conclude  this  section  and  the  article  with  the  remark 
that  the  relation  (56)  for  V  reduces  to  a  simple  expression  in  the  case 
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in  which  V  is  represented  by  the  kernel  v(x,x'  )•   This  expression  gives 
v(x,x'  )  directly  in  terms  of  q(x,x' )  which  is  a  solution  of  the  integral 
equation  (71).   Fran  (56)  we  obtain 

V  =  (AK-  KA)K  , 
or  equivalently 

v(x,x'  )  =  A{q(x,x'  hW  -   x)}- A'  (q(x,x'  )n(x'-  x)} 

+  /  [A{q(x,x")Ti(x"-x))]  {q^(x",x)il(x- x")}dx" 
-  j^A"[{q(x,x")Tl(x"-x)}](q^(x",x)Ti(x-  x")}da". 

If  we  express  the  functions  and  Laplacians  on  the  right  in  spherical 
coordinates  then  because  of  the  fact  that  the  derivative  of  the  t^ 
function  is  a  &  function  we  find  that  v(x,x'  )  is  a  smi   of  two  terms,   ' 
one  of  which  has  &(x'-x)  as  a  factor  and  the  other  of  which  has  t)(x'-x) 
as  a  factor  (the  terms  in  6'(x'-x)  cancel  out).   Since  V  must  be 
Hermit ian  we  have 

V  (x' >x)  =  v(x,x'  ), 

and  thus  the  coefficient  of  the  term  in  ti(x'-x)  must  vanish.   What 
remains  is  the  following  expression  for  v(x,x'  ): 
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(72)  v(x,x')=  ^ ^^. 


2 
X  q(x,0,cp;  X, 


where  q(x,0jCp;  x,9' ,cp'  )  is  obtained  from  q(xjX' )  by  using  spherical 
coordinates  and  setting  x' =  x.   We  recall  that  q(x,x' )  is  the  solution 
of  the  integral  eqiiation  (71)  whose  kernel  w(x,x'  )  is  the  representer 
of  the  operator  n  =  W-  I  (cf.  (53)  )  • 


APPENDIX  I 


Proof  of  the  Identity 


±ik|x-x'|    ,.     -,     p   i£- (x"x'  ) 
1  e    '-  ~  '  _  lira    1     /  e  -^ 

hn      |x-  x'  I    ~  ^  ^  °  (; 


/ 


3  J  ,2 


d£ 


We  have 


k  -  p  ±  ie 


i£.(x-x')        .    p  p  p  „ipU-x' [cos 


i-r  /  % 5 d£  =  -^  /  /  /  ^-2 2 P  "^^  ®  ^P^®^^ 

l^)^J    k^-p^iie      (2«)^^„^M    k^-p^±ie 


o  o  o 


r  p-ipU-x'  I  „ip|x-x'l 

^^ / Pdp 


(2«)^|x-x'ro      k^-p^±ie 


p  sm  p 


|x-x'  I  1        r  P  sin  p|x-x'  I 


P  P  faJ-n  p|x-i  I 

dp  = ^; — r— ^ ^P 


(2rt)^|x-x'|  ^   k^-p^±ie        (2n)^|x-x'|  _^   k^-p^±ie 
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2i(2jt)    |x-x'  I 


P             iplx  -  x' I       , 
/  e  ^' 'p  dp 


/ 


•ip|x  -  x'  I 


P  (ip 


o    (/?±ie-pj(yi?±Te+p)    -co     {/k^±ie--p)[s/k^±  ie  +  p) 


1 

.^e±ik|x-x'l.  ig  ik|x-x'I 

1 

I|7t|x-X'| 

-"   ^n|x-x'| 

APPENDIX  2 

Proof  of  tl 

le  asymptotic  relation 

±ik|x  -  x' I 


ikx 


(A.l)     ^^^_^  e  -"Vsin  e  ~  -2ijTkS(e  -0'  )S(cp-  cp'  )  ^ 


+2i7tk6(jt  -  e  -  e'  )6(cp  ±  jr  -  cp'  )  — 


-ikx 


We  begin  with  a  proof  of  the  relation 


+  •      2 


±^2 


, .    ^  V        lim       2"     ±iax  n  "^   ■^t    \ 

{k.2.)     ^     _  CL''   e  =    x/jt  e  &(x). 


Proof  of   (A.2): 

Consider 

°°2  1^2  \      "  2. 

lim  2    /      ±ia-x   „/    >,  lim       2    T    ii-°^x    _/    ^,      ,      lim       2.    ?    ±iax    „/    ^^ 

^   ^        a''    /   e  f(x)dx  =  a''    /    e  f(x)dx   +  a         e  f(x)dx 


2    r    ±iax   „/    > ,  lim       i    f 

I.      /   e  ffxjdx  =  a'^    / 

o 

■       p      ±it  p     ±: 


f(0)   =  y^  e       ^  f(0), 
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where  we  have  used  the  change  of  variable 


The  proof  of  (A.l)  can  be  carried  out  by  the  method  of  stationary  phase. 
We  write  first 

k-x  =  kx[cos(cp  -cp'  )sin  0  sin  G'  +  cos  0  cos  0'  ]. 

Stationary  points  where  derivatives  of  k-x  with  respect  to  0  and  cp 
vanish  are 

0'  =  0,  cp-!  =9 


cp  +  jt,      0<cp<rt 
cp  -   jt,      n  <  cp  <  2rt 


We  consider  separately  two  power  series  expansions  of  k-x  up  to 
quadratic  terms  about  each  of  these  stationary  points  in  the  0'  ,cp' 
plane.   The  expression  e  -  "%  sin  0  can  be  regarded  as  a  sum  of  two 
terms  in  which  the  phase  k-x  of  the  exponential  is  replaced  by  the  power 
series  about  each  of  the  stationary  points  in  turn.   If  this  is  done 
and  (A. 2)  is  used  we  arrive  at 
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lim   ik-x,2  .   _    lim   ikx 
e  -  "^  sin  0  =      e 
X  ->  »  X  ->  -» 


^  ^   /xk  sin  9  ^  ^/xk  sin 


X  k  sin  9  5(9-9'  )S(cp  -  cp  ) 


lim   -ikx 
e 


^-   -  'yri^^  =  y^dl^  '^'''"   e'B(,-9  -8')6«p±„  V) 


ikx  g-ikx 

-2ijtk6(9-0')&(cp  -cp'  )  ^-—  +   2i«k&(jt-9  -9'  )&(cp±rt-cp'  )  —^ 
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